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The theoretical description of interacting fermions in one 
spatial dimension is simplified by the fact that the low en- 
ergy excitations can be described in terms of bosonic degrees 
of freedom. This fermion-boson transmutation (FBT) which 
lies at the heart of the Luttinger liquid concept is presented 
in a way which does not require a knowledge of quantum field 
theoretical methods. As the basic facts can already be intro- 
duced for noninteracting fermions they are mainly discussed. 
As an application we use the FBT to present exact results 
for the low temperature thermodynamics and the occupation 
numbers in the microcanonical and the canonical ensemble. 
They are compared with the standard grand canonical results. 



I. INTRODUCTION 

The low temperature thermodynamic properties of 
simple metals can be qualitatively understood in terms of 
the simple Sommerfeld modelcl, which treats the conduc- 
tion electrons as noninteracting fermions in a box. Typi- 
cal results are a specific heat linear in temperature and a 
constant spin susceptibility in qualitative agreement with 
experiments. Landau's Fermi liquid theoryn rests on the 
assumption of quasi-particles which are in a one-to-one 
correspondence to noninteracting fermions. This leads 
to a linear specific heat and a constant spin susceptibil- 
ity but involves renormalized quantities like the effective 
mass and the quasi-particle interaction parametersB'0, 
which arc difficult to calculate microscopically. The con- 
sistency of the approach was shown using perturbation 
theory to infinite order and more recently by renormal- 
ization group techniquesu. 

The problem of interacting fermions simplifies in one 
dimension. In a pioneering papem Tomonaga treated the 
case of a two-body interaction which is long ranged in real 
space. He showed that the low energy excitations of the 
noninteracting as well as the interacting system can be 
described in terms of noninteracting bosonsB. The im- 
portant idea to solve the case of interacting fermions was 
the observation that a long range interaction in real space 
is short ranged in momentum space and therefore only 
particles and holes in the vicinity of the Fermi points are 
involved in the interacting ground state and states with 
a low excitation energy. To obtain his results, Tomonaga 
linearized the energy dispersion around the two Fermi 
points ztkp. Luttingertl later used a model with strictly 
linear energy dispersions_which is closely related to the 
massless Thirring modelo. The exact solution for the 



Luttinger model was presented by Mattis and Liebtl. A 
very elegant method to calculate correlation functions 
for the pncLel is the bosonization of the fermion field 
operatoroO. The exponents of the anomalous power 
law decay of various correlation functions are determined 
by the anomalous dimension, which can be calcula 
explicitly for the Tomonaga-Luttinger (TL-) ppdeflt 
It was an important observation of HaldaneO that the 
low energy physics of the exactly solvable TL-model pro- 
vides the generic scenario for one-dimensional fermions 
with repulsive, interactions. Like in the Landau Fermi 
liquid pictureel a few parameters completely determine 
the low energy physics. Generally they are as difficult 
to calculate as the Landau parameters. In contrast to 
the higher dimensional case there are additional exactly 
solvable models for which Haldane's Luttinger liquid sce- 
nario can be tested and the parameters determining the 
anomalous dimension can be calculated using the Bethe- 
ansatz techniquetHl. 

Another important manifestation of Luttinger liquids 
is called spin charge separation, i.e. for low energy ex- 
citations the charge and spin degrees are completely de- 
coupled. This shows up, for example, in the .spectral 
function of the one-particle Green's functionE^O which 
largely determines the photoemission spectrum. Recent 
high resolution photoemission experiments on quasi one- 
dimensional organic conductors have been interpreted to 
show Luttinger liquid behavioura. 

In the present paper we present the basic ideas be- 
hind the Luttinger liquid concept in very simple terms, 
involving only basic quantum mechanics and statistical 
mechanics. We therefore believe that our approach re- 
quires less theoretical knowledge than the field theory 
appxoach to FBT published in this journal some time 
agdij. 

In section II we present the fermion-boson transmuta- 
tion in terms that should be accessible to anybody who 
has had an elementary course in solid state physics. Es- 
pecially no use is made of the method of second quanti- 
zation. The concepts introduced in section II are used in 
section III and IV to compare the different statistical en- 
sembles for noninteracting fermions in a boxliZHij. In the 
grand canonical ensemble the constraint of a fixed parti- 
cle number which makes the calculation complicated in 
the canonical ensemble is relaxed and the determination 
of thermodynamical quantities is straightforward. The 
microcanonical ensemble is seldom used in this context. 
Our discussion of the basic ideas of the FBT allows us to 
present a comparison of the three ensembles. We show 



1 



explicitly how the results approach each other in the limit 
of infinite box size with constant density, i.e. the ther- 
modynamic limit. For the calculation of the occupation 
numbers in the canonical ensemble the numerical effort 
increases quickly with the system size when the elemen- 
tary approach of section IV is used. We take this as 
an opportunity to discuss the fermion-boson transmuta- 
tion again on a more advanced level using the method of 
second quantization in section V. The techniques which 
allow the calculation of-,the properties of the TL-model 
of interacting fermionsou are introduced again in the 
context of noninteracting fermions. The power of this 
method is shown by presenting a closed expression for 
the occupation number of noninteracting fermions in the 
canonical ensemble, which allows to discuss the transition 
to the Fermi function in the thermodynamic limit. 

The paper is written such that readers not familiar 
with the method of second quantization can understand 
fcrmion boson transmutation in its elementary form and 
the comparison of the statistical ensembles by reading 
sections I to IV only. 

II. NONINTERACTING FERMIONS IN A BOX 

We consider particles in one dimension in a box of 
length L with fixed boundary conditions 0(0) = 4>{L) = 
0. The single particle eigenfunctions and energies are 
given by 

4>n(x) = \/2/Lsm(k„x), (1) 

and 

e n = e(k n ) = -^, (2) 
2m 

with fc„ = nn/L and n G IN. Later we will also con- 
sider periodic boundary conditions. Then the spacing of 
the A;-values is doubled and negative and positive values 
occur. In the following the spin of the fermions is ne- 
glected for simplicity, i.e. we consider spinless fermions 
as in [ HQ. The ground state of N fermions is given 
by filling the lowest momentum states up to the Fermi 
wavevector kp — npir/L where np = N and the ex- 
cited states can be classified by giving the occupation 
numbers Hk n , which are zero or one. For fixed particle 
number N the sum of the occupation numbers has to be 
equal to N. This constraint makes the evaluation of the 
canonical partition function difficult. In the well known 
procedure of Darwin and FowlerES an integral represen- 
tation for the Kronecker delta describing the constraint 
is used. If the corresponding integral is evaluated by the 
method of steepest descent the result is equivalent to the 
grand canonical calculation, where only the average par- 
ticle number is fixed as particle exchange with a bath is 
assumed. In the following we show that in the low tem- 
perature regime the canonical partition function can be 



expressed as the partition function of a photon gas, i.e. 
massless bosons. 

We .first summarize the well known grand canonical 
resultsED. The total energy E gc as a function of the tem- 
perature T and the chemical potential fi is given by 

n 

where f3 = l/(fcsT) and fcs is the Boltzmann constant. 
The chemical potential /i becomes a function of tempera- 
ture if the average particle number is assumed to be given 
by the fixed value N 

N = H e /3[ £ „- M (T)] + 1 - ( 4 ) 

n 

In the thermodynamic limit L — > oo, N/L fixed, the low 
temperature result for the energy, the specific heat and 
other thermodynamic quantities can be-explicitly evalu- 
ated using the Sommerfeld techniqueElO. They are de- 
termined by the one-particle density of states per unit 
length po(e) = l/i:hv(e) in the vicinity of the Fermi en- 
ergy ep = e np , where % times the velocity v(e) is given 
by de(k)/dk as a function of e. The low temperature 
specific heat, for example, is linear in temperature and 
determined by p(e F )tEl This is an indication that in 
order to arrive at this result one is allowed to simplify 
the spectrum of eigenvalues in Eq. (2) in the vicinity of 
the Fermi energy. 

For n 3> 1 the eigenvalue spectrum given by Eq. (2) 
becomes locally equidistant 
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FIG. 1. Specific heat for the quadratic (solid line) and the 
linearized (dashed line) electron dispersions. 

As in Refs. [ § and [ § we therefore linearize the 
spectrum e„ around the Fermi point np 
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ej/ = ep + hvp{n — uf^/L, (6) 

where vp = hkp/m = Kiz(N/L)/m is the Fermi velocity. 
As we keep the density (N/L) constant while increasing 
L the Fermi velocity is independent of the system size. In 
the following we consider this as the model for which we 
want to determine the low temperature thermodynamics. 
Fig. 1 shows the temperature regime fc^T <tf for which 
the specific heats for both models agree. The results 
where obtained by numerical evaluation of Eqs. (3) and 
(4) in the thermodynamic limit. The linearization in Eq. 
(6) is the essential step ta-obtain the exact solution for 
the interacting TL-modcloLl The constant energy sepa- 
ration of the linearized model is given by A = hvpir/L. 

For the linearized model also the excitation spectrum 
of the TV-electron system is equidistant 

E M = E + MA, (7) 

where Eq is the ground state energy of the filled Fermi 
sea. For large M the eigenvalues are highly degenerate. 
The determination of the degeneracy can be reduced to 
an old problem in combinatoricstJ, if we classify the cor- 
responding states not by fermionic occupation numbers 
as discussed above, but in terms of the upward shifts of 
the occupied levels compared to the Fermi secB. This is 
exemplified in Fig. 2 where a special excited state for 
M = 20 is shown. The highest occupied level is shifted 
by seven energy units, the second and third by four, the 
forth by three and the fifth and sixth highest by one unit. 
All deeper levels are unshifted. The excitation energy in 
units of A is given by M = 7 + 4 +.4-^3 + 1 + 1 = 20. 
This is one of the possible partitionstnEl of the number 
M = 20. 
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FIG. 2. Example of the classification scheme for the excited 
states. 

Each other of the 627 partitions of M = 20, written in 



lexicographic order, as in our example, uniquely classi- 
fies an other possible excited state corresponding to this 
excitation energy. For small M it is easy to write down 
all partitions. For M = 4 the five possible partitions are: 
[4]; [3, 1]; [2, 2]; [2, 1, 1]; [1, 1, 1, 1]. For a general M a par- 
tition is defined by a set of M numbers m, € INo such 
that 

M 

M = mi + 2m 2 + 3m 3 + ... + Mm M = ^ Imi (8) 

i=i 

Already more than two hundred years ago EulerEH proved 
important theorems about the number of partitions. His 
use of generating functions can be considered the earliest 
precursor of the methods of quantum statistical mechan- 
ics. 

A formally simple expression for the number gM of 
partitions of M can be given as a sum over a Kronecker 
delta using the defining Eq. (8) 

oo oo oo 

mi— niM=0mM+l=0 

From Eq. (8) it is obvious that it is unnecessary to include 
more than M summations in Eq. (9), but the form with 
the additional infinite number of summations turns out 
to be useful to calculate the canonical partition function 
of the linearized model 

oo 

Z C =J2~9^ (E ° +MA) - ( 10 ) 

M=0 

Here gM is the degeneracy corresponding to the excita- 
tion energy MA. For excitation energies MA smaller 
than up A all partitions of M lead to a possible ex- 
cited state of the /V-particle system and gM = 9 hi- For 
M > up the degeneracy of excited states is smaller than 
the number of partitions ()m < 9 hi- For M = np + 1, 
e.g. there is no allowed excited iV-particle state related 
to the partition 

[1,1, 

np+l 

because in the groundstate there are only uf fermions 
that can be shifted upwards. For M > np we have to 
take into account that one-particle states exist only for 
quantum numbers n > 1. Two arguments can be given 
to use gM as the degeneracy also for M > np. The one 
in spirit of Tomonaga'scl treatment of the interacting case 
is to argue that for n_pA/(fcsT) 3> 1 the contribution of 
these highly excited states is exponentially small in this 
large ratio and the replacement of c/m by 9m in Eq. (10) 
makes no relevant difference. The solution corresponding 
to Luttinger'sQ treatment of interacting fermions consists 
of adding an infinite Dirac sea of fermions in one-particle 
states with n < 0. Then the replacement of cjm by gM 
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is exact, but one has to be careful with the treatment of 
an infinite number of electronsEftl 

After the replacement cjm — * 9m the calculation of the 
partition function Z c in Eq. (10) is simple if we use Eq. 
(9) 



oo oo 



-PE, 



xe 



E E - E 

M=0mi=0 m M =( 



Vy° lr, 



n e< 

1 = 1 Vrnj^O 
oo 1 



-miPlA 



(11) 



1=1 



This is the partition function of a system of massless 
bosons or harmonic oscillators with energy I A — hvpki. 
The nil which describe the partitions play the roles of 
bosonic quantum numbers. This is FBT in its most ele- 
mentary form. 

If we differentiate — In Z c with respect to f3 we obtain 
the canonical energy E c 



E C (T) = E + J2 



I A 



(12) 



The evaluation of Eq. (12) for finite systems will be 
discussed in the next section. Here we first want to 
show that the results for the energy per unit length 
e(T) = [E(T) - E ]/L using Eq. (3) and Eq. (12) agree 
in the thermodynamic limit. In this limit the sums can 
be replaced by integrals and we obtain from Eq. (12) 



e c (T) = — r 1-, 

Trnvp 

where the integrals I± are defined as 



/+ = 



e x ±l 



-dx. 



(13) 



(14) 



Eq. (13) is the-<pne-dimensional version of the Stefan- 
Boltzmann lawcll of black body radiation which states 
that in d dimensions the energy density is proportional 
to T d+1 , or to Debye's theory of theJow temperature 
specific heat due to acoustical phononstl. 

For the calculation in the grand canonical ensemble us- 
ing Eq. (3) the contribution of the states above ef equals 
the part below tp for the linear dispersion and the chem- 
ical potential is temperature independent. This yields 



{k B T) 2 
■kTivf 



2J 4 



(15) 



The integrals I± can be found in every good table of 
integrals or textbook on statistical mechanicsEZI. The fact 
that the results in Eqs. (13) and (15) agree can be shown 



without actually calculating the integrals, 
algebra yields 



Elementary 



I--I+ 



-dx = 7_/2. 



(16) 



and we obtain J_ = 2/ + , which shows the equivalence of 
the canonical and the grand canonical ensemble in the. 
thermodynamic limit. The numerical value is given bycll 



C(2) 



"/6 where Q{x) is the Riemann zeta func- 



tion. The specific heat per unit length cl = Cl/L for 
both ensembles is therefore in the thermodynamic limit 
given by 



irk 2 B T 

3 flVF 



(17) 



In the following section the specific heat for both ensem- 
bles is discussed for finite systems. For a given density 
and temperature there is a natural length scale (thermal 
de Broglie wave length) given by a(T) = Tivfk/ (hsT). 
The system size L is large compared to a(T) if k B T is 
large compared to the level spacing A. 

So far we have not discussed the microcanonipal en- 
semble in which one works with a fixed energyO and 
the entropy as a function of energy is determined by the 
logarithm of the degeneracy. Therefore we have for the 
linearized model 



Sn 



,(E M ) = k B In (sm), 



(18) 



where gu is the number of partitions. There are recur- 
sive methods to determine <?m for small and intermediate 
vahies-pf M and excellent estimates can be given for large 
Afc3c3. We show also in the following section how the 
microcanonical results approach those of the canonical 
and the grand canonical ensemble for large L. 



III. THERMODYNAMIC PROPERTIES FOR 
FINITE SYSTEMS 

The Stefan-Boltzmann law Eq. (13) for the energy per 
unit length is only correct in the thermodynamic limit. 
For finite L the sum in Eq. (12) has to be performed nu- 
merically. For L/a(T) ^.1 the Eulcr summation formula 
in its most simple forrr£3 

E /(?') - r f( x ) dx + [/Oo) + /(ii)]/2 T (19) 

7=1n J 3« 



provides the leading finite size correction 



Cl,c = 1 3 a(T) 

Coo 27T 2 L 



a(T) 



L 



(20) 



Fig. 3 shows cl, c /coo as a function of L/a(T). As 
L/a(T) = ksT/ A this plot can be viewed as the specific 
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heat ratio as a function of system size for fixed tempera- 
ture or as a function of temperature for a given length of 
the system. For small arguments the ratio is small due 
to the well known phenomenon of freezing of the excited 
states far_temperatures small compared to the excitation 
energiesO. Fig. 3 also shows the corresponding results for 
the grand canonical specific heat following from the nu- 
merical evaluation of the result using Eq. (3). For small 
arguments the exponential suppression is weaker than in 
the canonical case as the chemical potential lies in the 
middle between the highest occupied and the lowest un- 
occupied state of the Fermi sea. The use of the Eulcr 
summation formula Eq. (19) shows that there is no finite 
size correction of order a(T)/L as in the canonical result 
Eq. (20) and the asymptotic limit is reached very quickly. 
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FIG. 3. Ratio of the grand canonical specific heat and the 
specific heat in the thermodynamic limit Coa (solid line) and 
the canonical specific heat and (dashed line) as a function 
of L/a(T). 
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dE 



A I 



^[ln 5 M]'. (22) 



Due to the discreteness of the energy the derivative 
in Eq. (22) has to be evaluated as a ratio of differ- 
ences f' M = (/m+i — ,/m-i)/2- This definition produces 
smoother results than the definition involving differences 
of neighbouring numbers. The microcanonical definition 
of the heat capacity for constant L is 
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L,micro 
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(E M ) 



dT„ 



dE 



M 



(23) 
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FIG. 4. Microcanonical entropy as a function of the energy 
calculated with help of the exact degeneracies (circles) and 
the Ramanujan approximation Eq. (18) (solid line). 
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As discussed in the preceding section the number of 
partitions gu plays a central role for the thermodynam- 
ics of the linearized model. In the microcanonical en- 
semble the number enters explicitly in Eq. (18). The 
values of the gM fp; M up to 500 can be found in a 
standard handbookej. They increase very quickly with 
M as can be seen from 3500 ~ 2.3 • 10 21 , which leads to 
S mi cro{E + 500A)/fc B w 49.20. For values of M of the 
order fifty q r ^ aij ger already the crudest approximation of 
Ramanujar 



(9m)r = e"V™t*/(4My/3), 



(21) 



provides a good approximation. One obtains, for exam- 
ple ln{g M (500)}^, w 49.18. In Fig. 4 we show S mlcro 
as a function of M and compare the exact result using 
the recursion relatiorcj'E 2 ] with the Ramanujan approxi- 
mation Eq. (21), which smoothes the number theoretical 
fluctuations. The inverse microcanonical temperature is 
defined as 
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FIG. 5. Ratio of the exact microcanonical specific heat and 
the one calculated with the help of the Ramanujan approxi- 
mation as a function of [(Em — Eo)/ A} 1 ^ 2 . 

For M>lan analytical approximation for this quan- 
tity can be obtained using the Ramanujan expression Eq. 
(21), which we denote as Cl.r 
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(24) 



This should be compared with the relation Lc^ — 
Trk By /2Le(T)/(3A), where e(T) = e c (T) = e gc (T), 
which follows from Eqs. (13) or (15). As in the large 
energy limit E — E in the microcanonical ensemble and 
E{T) — Eq in the canonical (or grand canonical) ensem- 
ble can be identified, as discussed in the next section, 
Cl.r/L agrees with c JOO in this limit. Fig. 5 shows how 
the exact microcanonical result Eq. (23) approaches Cl,r 
as a function of \J {Em — £^j/A, which takes the role of 
L/a(T) in Fig. 3. The ratio 

cl, micro I ' c l,r shows quite 
strong fluctuations for values yj (Em ~ Eq)/ A < 5. 

As in the thermodynamic limit the entropy per unit 
length agrees for all three ensembles we can use the grand 
canonical entropy to obtain an estimate for the number 
of partitions. This yields the Ramanujan expression Eq. 
(21) with 4Ma/3 replaced by 1. 



IV. OCCUPATION NUMBERS 

For variable particle number the method of second 
quantization in Fock space is the appropriate mathe- 
matical framework. The Hamiltonian for noninteracting 
fermions considered in section II reads 



„t 



(25) 



where the the occupation number operators and 

c m ( c m) are fermionic creation (annihilation) operators, 
which obey canonical anticommutation relations 



[c m ,C;]+ = 0, [c 



,tl 



(26) 



Explicit use of the method of second quantization is only 
made in the next section. Here we discuss the calcula- 
tion of the expectation value of n m in the three ensem- 
bles. The-calculation in the grand canonical ensemble is 
standardO and for arbitrary system size yields the Fermi 
function 



1 



Men 



-M) + 1 



= /(em). 



(27) 



We next consider the model with the linearized disper- 
sion Eq. (6) in the microcanonical ensemble. In section 
II we have labeled the states of excitation energy MA in 
terms of the list of upward shifts [a] m 



[a\M 



h,h, ■••] with h > h+i > 0, 



(28) 



which corresponds to a partition in lexicographic order. 
The expectation value of the occupation number is given 
by 



' micro, h/M 



— Y] ([a}Ai\n m \[a} M ) 



9m 



(29) 



_a] M 



The expectation value of n rn in the state | [o]m) can easily 
be calculated due to the definition of the list [cy\m 

(30) 



1, 



( [a] m | rim | Mm) = 1 for m = n F + h - 

and zero otherwise. In order to perform the numerical 
calculations all partitions of M have to be explicitly cre- 
ated, which leads to a drastic increase of computer time 
for increasing M. The numerical results will be discussed 
below together with the canonical result 



b ra/ c 



Y gM (™ m ) 

c M=0 



ro,E h 



-0(E o +MA) 



(31) 



The canonical result can therefore be obtained directly 
from the microcanonical one if the temperature is low 
enough that the probability distribution 

w M (T)=g M e-K E °+ MA )/Z c , (32) 

has decayed to "zero" for values of M for which the calcu- 
lation of the microcanonical average is still feasible. Note 
that while the Fermi function can be obtained from the 
graad canonical potential by a derivative with respect to 
e m EZI such a procedure is not possible for the canonical 
free energy. In order to calculate the canonical average 
for higher temperatures, the method of the bosonization 
of the field operatorEjO can be used which is described 
in the next section. We therefore only consider results 
for quite low temperatures in this section. 

Fig. 6 shows the average occupation numbers for the 
three ensembles for r = hst/A = 3. The microcanonical 
curve is for M — 10, which roughly corresponds to the 
maximum of the probability distribution wm(T) shown 
in Fig. 7. The overall agreement of the results is quite 
good already at this small value of r. Note that the 
microcanonical result yields the ground state occupancies 

for \m-n F \ > M, e.g. (n m ) micro . EM = for m > M + n F, 
while the canonical and the grand canonical results go to 
zero continuously. 
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FIG. 6. Occupation numbers for the three ensembles for 
r = 3 respectively M — 10. 



5 



0.04 



0.02 



0.00 




FIG. 7. Probability distribution of the energy for 
t = k B T/A = 3. 



V. BOSONIZATION OF FERMIONIC 
OPERATORS 

This section is more technical than the previous ones 
and some familiarity with the method of second quanti- 
zation is required. In section II we have shown that the 
N fermion states of the model with the linearized dis- 
persion are highly degenerate and can be labeled by the 
list of upward shifts corresponding to partitions, which 
leads to the bosonic canonical partition function. We now 
show that one can choose linear combinations of the de- 
generate states for a given M, which can be described by 
bosonic creation and annihilation operators (ladder oper- 
ators) . They play the essentiaLeole for the solution of the 
interacting fermion problemtlO'El. One defines operators 
bilinear in the fermionic operators 



1 

Vl 



71=1 



(33) 



with I > 1. This is a sum of operators which all shift a 
fermion upwards by I energy units. Using the fermionic 
anticommutation rules it is easy to show that [b,, b,,] =0 
and 



b f ] - - V J c 



(34) 



The operator on the rhs counts the number of fermions 
in the lowest I states at the bottom of the band divided 
by I. For I <C tif and N fermion states which have no 
holes at the bottom of the Fermi see one can replace the 
operator on the rhs of Eq. (34) by the unit operator. One 



can make this an exact statement if a Dirac sea is added, 
i.e. the sum in Eq. (33) runs not from 1 but from —mo 
with mo — > oo. An analogous discussion for I ^ V leads 
to 



[bi 



4] = w 



(35) 



i.e. the operators defined in Eq. (33) obey boson commu- 
tation relations. 

With the help of these ladder operators we can create 
an orthonormal basis set of many fermion states labeled 
by bosonic quantum numbers {m,;} with rhi € INo, if one 
treats the Fermi sea \F(N)} as the harmonic oscillator 
ground state 



iw)=n- 



i 



(b\r \f(n)) . 



(36) 



All states with Irhi — M provide a new basis of states 
corresponding to the excitation energy MA of Hq. They 
are rather complicated linear combinations of the states 
discussed in section II which could be labeled either by 
fermionic occupation numbers or by the list [o\m of the 
upward shifts. If for M = 4 we label the |[o:]m) by \i) 
with i = 1, ..,5, in the order in which the partitions of 
M = 4 were presented in section II, we obtain e.g. 

b\\F) =i(|i)-|2) + |4)-|5)) 



b\b\ \F) 



V4! 



(b\) 4 \F) = 



1 
1 



(|1>-|3> + |5» 
(|1)+3|2) 



+2|3)+3|4} + |5». 



(37) 



As long as we consider only noninteracting fermions both 
types of basis sets are completely well suited for the de- 
scription. It turns out that for fermions with a two-body 
interaction the interacting groundstate has a much sim- 
pler form in the bosonic basis set Eq. (36). 

In Eq. (33) we have expressed the Bose operators h{b\ ) 
in terms of the fermionic operators. In order to take ad- 
vantage of the bosonization one has to address the reverse 
problem: Can one express fermionic operators completely 
in terms of the Bose operators'? Here one proceeds in two 
steps. It is possible to bosotuzc the operator Hq with the 
help of the Kronig identityliil 



i=i 



E 

n=l 



ncic n --{N 2 + N), 



(38) 



where N = ^2 n c^c„ is the fermionic particle number 
operator. For the proof of this relation the boson com- 
mutation relations of the bi(b\) is not needed. As the lhs 
of Eq. (38) contains products of four fermion operators 
it might look surprising that one obtains a one-particle 
operator and only a term involving the particle number 
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operator on the right hand side. It is left as an exercise 
to the reader to show that the additional four fermion 
terms exactly vanish. If we subtract npN and multiply 
by HvFTt/L we obtain 



L 

1=1 



b}bi + h(N), 



(39) 



where the function h(N) is irrelevant in the following. 

It as also possible to bosonize a long range two-body 
interaction if periodic boundary conditions are used. 
Generalizing the procedure described above one linearizes 
around both Fermi points ±fcp. Then the Bose operators 
corresponding to Eq. (33) are the two contributions to 
the Fourier transformed density operators. As the two- 
body interaction is quadratic in the density the two-body 
operator is bilinear in the Bose operators. The problem 
of interacting fermions is therefore reduced to a quadratic 
form in Bose operators, which is easily diagonalized by 
a canonical transformationfl'EI We do not further dis- 
cuss interacting fermions here but want to calculate the 
expectation value of n m in the canonical ensemble for 
noninteracting electrons. For that purpose we have to 
express n rn in terms of Bose operators. In order to do 
that we have to proceed similarly as in the bosonization 
of the field operatoi^tli. 

This requires some preliminary steps. Let us consider 
the following operators linear in the boson operators 



E< 

n=l 



k„&+ ; B_ = (3 n b n , 



(40) 



where a n and (3 n are arbitrary constants. Using the bo- 
son commutation relations Eq. (35) one can prove the 
commutation relations 



[b m ,e 



G 1 — G-jYi G G 

[bl,e B -e A +]=-p m e B -e A < 



(41) 



If we can construct an operator S in terms of the 
fermionic operators which fulfills the commutation rela- 
tions [b m ,S] = —a m S and \b' m , S] = f3 m S then Se B ~e A+ 
and 6^ and can therefore only be a 
operator N, which commutes 



commutes with all b r , 



function of the numb 
with the boson operators^, i.e. 



S 



G N e ~ 



(42) 



In order to bosonize c}jC n we consider the more general 
operator etc; with £ Z , i.e. the infinite Dirac sea is 
included, and determine its commutation relations with 
b m and b] n . Using the anticommutation relations Eq. 
(26) one obtains 



[bm,c{ci] = _( c t_ m q - c f k c l+m ) 

-7=( c l+ m c l - 4 c *-™)- 



(43) 



If we now multiply Eq. (43) by e ~ lku e llv with real u and v 
and sum over k and I the commutation relations with the 
resulting operator are of the form as discussed following 
Eq. (41) 

[b m ,^(u)i>(v)] = -±={e- imu - e- imv )^(u)i>(v) 



[bl,tf(u)1>(v)] = —(e imu -e imv )^(uMv), (44) 
and the auxiliary field operator ip(v) is defined as 

oo 



l=-c 



It is periodic in v with period 2-7T. Using the argu- 
ments presented following Eq. (41) the bosonized form 
of ^(u)^p(v) reads 



where the operator <f>(u) is given by 



(u) = -iV^n, 
n— 1 v 



(46) 



(47) 



and Gjy(u, v) is a so far undetermined function of the par- 
ticle number operator. As the Fermi sea with N fermions 
in addition to the Dirac sea is the vacuum state for the 
bosons, i.e. 6; \F(N)) = 0, the expectation value of the 
product of the exponentials in Eq. (46) in the Fermi sea 
equals one and yields 



G N (u,v) = (F(N)\^(u)i>(v)\F(N)) 



(48) 



which can be calculated in the fermionic picture using the 
definition Eq. (45) and ^F(A)|cj^c ; \F(N)\ = S kt iQ(N - 
I), where Q(x) is the step function (8(0) = 1) 



G N {u,v) = ^2 



Jl(u-v) 



/ — — oo 



E 



a il(u— v— iO) 



(49) 



In the second equality we have added an infinitesimal 
imaginary shift in u — v in order to make the formal sum 
convergent. This is a typical example of the problems 
with the infinite Dirac sea. 

After the bosonization of ip'(u)ip(v) is completed, 
which actually is simpler than the bosonization of a sin- 
gle field operator as no additional operators changing the 
fermion number have to be introducecO, we can come to 
our application of calculating the canonical expectation 
value of c^c m . If we denote the canonical average of an 
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operator A in a state with N fermions in addition to the 
Dirac sea by (A) NT we obtain with the Bose function 
b(x) = (e x - l)- 1 



(u)ip(v)) N ^ T = G N {u,v) 
x cxp | ( 2 cos i n ( u - v )] - 2 ) b ( nA ^ 



(50) 



Here we have used the formula (e A e B ) = 6 H a2 + 2AB + b2 ) 
for a canonical expectation value for free bosons which 
was proved by MermirEJ in "one sentence" . It is valid for 
operators A and B which are linear in the boson creation 
and annihilation operators. To calculate the expectation 
value of the occupancies we need the inversion formula 
to Eq. (45) 



1 

2^ 



-i vl 



tp(v)dv. 



(51) 



As the rhs of Eq. (50) is a function of u — v only one 
integration can be performed and one obtains 



1 



c N+m C N+rn ) — 7T~ / e 

1 N,T 2lT J 



2tt 



—iu{m+N) 



x (V + («)V(0)) 



N.T ' 



(52) 



In order to analytically perform the remaining integra- 
tion in Eq. (52) we expand (ip^(u)if>(Q)) in a Laurent 
series in z = exp (iu) 



JuN 



N,T 



(53) 



n— — oo 



where the coefficients d n can— be calculated recursively 
as discussed in the appendixEil. Using the notation of 
section IV we finally obtain 



(54) 



For all values of r for which the calculation of (n m ) c T 
using Eq. (31) is numerically feasible (e.g. r = 3 in Fig. 
6), the results using Eqs. (31) and (54) are identical. As 
mentioned in the introduction the result of the present 
section which used Eq. (46) as the essential ingredient is 
very close in spirit to the calculation of the occupancies 
of the interacting TL-model. 
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FIG. 8. The difference between the canonical occupation 
numbers and the grand canonical Fermi function at r = 10 
(solid line) and r = 15 (dashed line). 

Fig. 8 shows the difference of the canonical and the 
grand canonical result for r = 10 and 15 as a function 
of m — uf ■ From a plot of the maximum deviation as a 
function of 1/t we find that it scales as 1/r, i.e. 1/L if 
the temperature T is kept constant. The relative devia- 
tion of the canonical occupation numbers from the Fermi 
function \[(n m ) c — f(e m )] / f{e m )\ for m > np increases 
with m — np and saturates in the exponential tail of the 
Fermi function. This shows that the canonical occupa- 
tion numbers decay as (1 — rj T ) exp { — (m — np — 0.5)/r)} 
for {jn — np)jT ^> 1. For r = 15, for example, r\ T = 0.032 
and rj T goes to zero like 1/r. In the thermodynamic limit 
L — > oo the sum in the exponent on the rhs of Eq. (50) 
goes over to an integral which can be performed analyt- 
ically. Performing the following integration in Eq. (52) 
one obtains the Fermi function, i.e. as expected and seen 
in Fig. 8 the canonical occupation numbers approach the 
grand canonical ones for L — > oo. 



SUMMARY 

This paper intended to serve a dual purpose. First 
to explain in terms as simple as possible the concept of 
bosonization which plays the essential role in the Lut- 
tinger liquid picture of interacting fermions in one di- 
mension. We only shortly discussed interacting electrons, 
but every reader who followed the algebra of section V 
should have no problem understanding e.g. the-,calcula- 
tion of the occupancies in the interacting caseEiJ. As an 
application of the concepts we presented results for the 
microcanonical and the canonical ensemble for noninter- 
acting fermions. The comparison of our results with the 
standard textbook grand canonical results can serve as a 
rare explicit example for the general arguments concern- 
ing the equivalence of the three ensembles in the thermo- 
dynamic limit. 
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APPENDIX A: 



In this appendix we present the method to obtain the 
coefficients d m of the Laurent series in Eq. (53) which de- 
termine the canonical occupation numbers (n nF+rn ) c T . 
According to Eqs. (49) and (50) we can write 



(^(u)m) N ,T = eiNu \J2 z ~ l j 

xexp{/(z) + /(l/z)-2/(l)}, (Al) 
where z — cxp (iu) and the function f(z) is given by 



n=l 



(A2) 



The geometric series Eq. (Al) is considered as a formal 
power series. The first step is to determine the coeffi- 
cients of the power series of the function 



F(z) = e'W = c mz m - 



(A3) 



m=0 



Using F'(z) = f'(z)F(z) and c = 1 one can obtain 
recursion relations for m > 1 



171 ' 



(A4) 



i=i 



The next step is to write the last factor in Eq. (Al) as a 
Laurent scries 



F(z)F{l/z)/[F{l)] 2 ^ £ ai z l . 



(A5) 



By comparing coefficients one obtains using Eqs. (A2) 
and (A3) 

{oo ^ oo 

-2 ^2 b( n / T )/ n \ X! c ™ c ™+i = a-i. (A6) 
n=0 J m=0 

As for fixed r the c m go to zero exponentially for large 
m the infinite sum converges well. In the last step we 
have to multiply the Laurent series in Eq. (A5) with the 
geometric series in Eq. (Al). This yields with Eqs. (53) 
and (54) 



(A7) 



71=0 



The numerical evaluation of the canonical occupation 
numbers with this method is possible for much higher 
temperatures than via Eq. (31). 
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